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Design details of a proposed high-energy (~50 kj/pulse), two-stage pulsed plasma thruster are pre- 
sented. The long-term goal of this project is to develop a high-power (~500 kW), high specific impulse 
(~7500 s), highly efficient (~50%), and mechanically simple thruster for use as primary propulsion in a 
high-power nuclear electric propulsion system. The proposed thruster (PRC-PPT1) utilizes a valveless, 
liquid lithium-fed thermal plasma injector (first stage) followed by a high-energy pulsed electromag- 
netic accelerator (second stage). A numerical circuit model coupled with one-dimensional current sheet 
dynamics, as well as a numerical MHD simulation, are used to qualitatively predict the thermal plasma 
injection and current sheet dynamics, as well as to estimate the projected performance of the thruster. 
A set of further modelling efforts, and the experimental testing of a prototype thruster, is suggested to 
determine the feasibility of demonstrating a full scale high-power thruster. 


Nomenclature 

A - cross-sectional area [m 2 ] 

C - thruster capacitance [F] 
hit - impulse bit [N-s/shot] 

I sp - specific impulse [s] 
j a - anode tip current density [A m~ 2 ] 
j QV - average current density [A in“ 2 ] 
k - Boltzmann constant [J K -1 1 
L - characteristic initial length [in] 

Lq - initial thruster inductance [H] 

m - mass [kg] 

mbit - mass bit [kg/shot] 

ni, n 2 - density [m -3 ] 

n e - electron number density [m 3 ] 
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1 Introduction 

Pulsed electromagnetic accelerators are devices which use 
intense bursts of electrical current ((9(10 4 — 10 6 ) A) to 
create high speed (0(1O 3 — 10 5 ) m/s) jets of plasma. They 
find application as plasma sources in many basic plasma 
science experiments[l] as well as in a specific genre of 
electric space propulsion device called the pulsed plasma 
thruster (PPT)[2]. The present work is motivated by the 
desire to improve the performance of pulsed electromag- 
netic accelerators in the context of plasma propulsion. 

In what follows, the design of a proposed high-energy 
PPT is presented. The emphasis of the presentation is con- 
ceptual; this paper serves to define the principles of oper- 
ation of the proposed thruster, provide rough calculations 
for the potential performance, and identify physical pro- 
cesses that will be the subject of future phenomenological 
modelling. We first provide a brief background of the state 
of PPT technology readiness - to serve as a backdrop for 
defining our motivation for pursuing the development of a 
new thruster. 
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1.1 Background 

Modern PPTs have the potential tor fulfilling the attitude 
control requirements on a satellite at greatly reduced mass 
and cost. They are also being considered for constellation 
maintenance for missions such as interferometric imaging 
of the Earth from space or deep space from an Earth orbit 
(c.f., Polzin et al[ 3]). The benefits of PPTs are their sim- 
plicity, very small impulse bits for precise control of satel- 
lite motion, reliability, and high specific impulse. Two 
classifications of PPTs exist, corresponding to the form 
of propellant used: gas-fed (GFPPT) or ablative propel- 
lant (APPT). The gas-fed variet) has the advantages of a 
“clean” exhaust plume and high specific impulse. The ab- 
lative version of the PPT uses a solid propellant, such as 
Teflon, to provide other advantages such as compactness 
and overall ease of system integration; however, plume 
contamination and lower specific impulse may limit the 
application of APPTs for some missions. 

From about 1960 to 1968 PPT research focused on the 
gas- fed variant (GFPPT). The GFPPT was envisioned as a 
potential “primary” propulsion system, where the GFPPT 
would process large amounts ol power (>100 kW) and 
provide enough average thrust (>1 N) to perform large 
At; maneuvers, such as interplanetary missions[4]. How- 
ever, by the late sixties, research turned toward steady- 
flow electromagnetic acceleration (the MPD thruster). 
The quasi-steady MPD thruster[5] allowed for the amorti- 
zation of efficiency-robbing losses intrinsic to pulsed op- 
eration (such as propellant loss at the leading and trail- 
ing edge of the gas injection pulse), through the use of 
a protracted current pulse and, hence, allowed for higher 
thrust efficiencies to be attained A “secondary” propul- 
sion niche, i.e., small At; attitude control maneuvers, still 
existed for the GFPPT. However, at about the same time, 
the ablative variety of PPT (APPT) was gaining favor, 
mainly due to its mechanical simplicity. The GFPPT re- 
quires the storage of gaseous propellant under high pres- 
sure and a fast-acting valve to meter puffs of propellant 
into the discharge chamber. Fuither, the propellant feed 
system in the GFPPT is required to operate in a leak- 
free manner for many (> 10 7 ) shots. A flight qualified 
system capable of performing these demanding tasks was 
not available at the time and, hence, only APPTs were 
flown[6]. 

This state of affairs remained until the mid-90’s, when 
the some of the negative issues related to propellant feed- 
ing in the GFPPT were ameliorated by the availability of 
reliable, fast valves and the development of high-speed 
solid-state electrical switching technology. The former 
development addressed the reliability issue, while the lat- 
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ter showed promise for substantially increasing the pro- 
pellant utilization efficiency. These technologies were im- 
plemented in the thrusters of Ziemer et al. [7, 8, 9]. These 
low energy (<10 J) thrusters show promise for replacing 
APPTs in missions which require small Av maneuvers. 
However, high-energy thrusters probably cannot benefit 
from the new switching technologies, since the necessar- 
ily high current levels preclude the use of solid-state de- 
vices. 

Even with the addition of these new technologies, PPTs 
remain one of the least efficient electric propulsion sys- 
tems (rjt <10%). However, the possibility of large incre- 
mental gains in performance make PPTs one of the most 
interesting devices, from a research perspective. 

1 .2 Motivation and Approach 

No one type of thruster is best suited for all types of mis- 
sions. For example, ion thrusters may be the best propul- 
sion option for a large Av , long duration mission; how- 
ever, if short trip times are required, their characteristi- 
cally low thrust density may preclude their use altogether. 
Similarly, Hall thrusters have attractive performance at in- 
termediate specific impulse levels (~2000 s), but funda- 
mental physical limitations prevent them from achieving 
the high specific impulse levels required for some mis- 
sions. In short, every type of electric propulsion device 
has a parameter space within which it performs best. It is 
our contention that the PPT may have unique capabilities 
to satisfy the propulsion needs for missions that require 
high specific impulse, high thrust efficiency, and high 
thrust density. The motivation of the present project is to 
develop a high power (~500 kW), high specific impulse 
(~7500 s), highly efficient (>50%) thruster for use as pri- 
mary propulsion in a high power nuclear electric propul- 
sion system. The thruster, which we call PRC-PPT1, uses 
a two-stage acceleration scheme, to circumvent some of 
the deficiencies which have plagued earlier thruster de- 
signs. 

High-energy PPTs have already been experimentally 
shown[10] to be capable of accelerating current sheets 
to speeds greater than 150 km/s (corresponding to a spe- 
cific impulse above 15,000 s, if 100% sweeping efficiency 
is assumed). Another study[4] demonstrated high-energy 
PPT operation at greater than 60% thrust efficiency (note: 
the efficiency was calculated using the estimated propel- 
lent mass inside the thruster during the current pulse; the 
cold-gas pulse was substantially longer in duration than 
the discharge and, hence, using the actual (total) cold-gas 
mass bit would have lead to a substantially lower calcu- 
lated value of thrust efficiency). These experiments prove 
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that there are no fundamental physical limitations that will 
prevent us from attaining our staled design goals. Let us 
state in one place, then, the practical limitations of exist- 
ing laboratory designs which must be addressed to bring 
high-energy PPTs to a level of maturity that will allow 
them to be considered as a viable near-term high-power 
electric propulsion option: 

1. Propellant utilization efficiency: Two mechanisms 
contribute to low propellant utilization efficiency in 
gas-fed PPTs. As described earlier, the delivery of 
discrete gas puffs to the discharge chamber using me- 
chanical valves is difficult. Ideally, a gas injection 
system would inject just enough gas to fill the dis- 
charge chamber and then quickly close. For typical 
thruster dimensions, a valve open-close duration of 
no more than about 100 ps is required to avoid pro- 
pellant leakage after the electrical discharge pulse. 
The design of a mechanical valve to meet this spec- 
ification is challenging. A second propellant loss 
mechanism in PPTs results from incomplete cou- 
pling of the Lorentz force density to the working gas. 
This inefficiency occurs when the accelerating cur- 
rent sheet is unable to fully entrain the propellant that 
it encroaches upon. Permeability of current sheets 
remains an active research topicfl 1]. 

2. Electrical efficiency: Typical PPTs suffer from poor 
electrical impedance matching between the source 
and the load. Oscillatory current waveforms result, 
which lead to undesirable "restrike” in the discharge 
chamber and significant energy dissipation inside the 
capacitor. 

3. Electrode erosion: Lifetime-limiting electrode ero- 
sion remains a major issue in assessing the via- 
bility of high-power PPTs. Electrode erosion has 
been more extensively studied in steady-flow plasma 
thrusters[12]. There remains a need for an experi- 
mental database of electrode erosion rates for PPTs, 
using a variety of electrode materials, propellants, 
electrode temperatures, cunent levels, etc.. 

4. Reliability: The reliability (and lifetime) of high- 
power PPTs is is impacted most strongly by the fol- 
lowing components: the gas injection valve, the ca- 
pacitors, the high current switch, and the (eroding) 
electrodes. New high-power PPT designs must ad- 
dress the reliability of each of these components, or 
somehow eliminate the need for them in the system. 

The "issues” that have been enumerated above define 
the challenges that confront the development of a new 


thruster, that is, new thruster designs should strive to pos- 
itively impact one or more of these PPT deficiencies. We 
now list, in a general sense, how our proposed thruster de- 
sign addresses each of these points (details of the design 
will be given in subsequent sections): 

L Propellant utilization efficiency: PRC-PPT1 uses a 
em liquid metal (lithium) feed system that requires 
no mechanical valve. Propellant loading is achieved 
through the vaporization and ionization of a liquid 
droplet (first stage); the size of the droplet is tailored 
to precisely deliver the desired propellant mass to the 
discharge chamber, where it is electromagnetically 
accelerated (second stage). Thus, problems associ- 
ated with using a gaseous propellant feed system are 
eliminated. To address the propellant sweeping ef- 
ficiency issue, PRC-PPT1 will operate in a higher 
gas density regime than is typical. Traditional PPTs 
have been plagued with poor efficiency because they 
have not been operated in a plasma regime that fully 
exploits the potential benefits of pulsed plasma ac- 
celeration by electromagnetic forces. PPTs have 
generally been used to accelerate low-density plas- 
mas. Operation of thrusters in this plasma regime 
allows for the development of certain undesirable 
particle-kinetic effects, such as Hall effect-induced 
current sheet canting[ 13]. PRC-PPT1 was designed 
to propel a highly collisional, dense plasma that has 
more fluid-like properties and, hence, is more effec- 
tively pushed by a magnetic field. Traditional PPTs 
also suffer from dynamic efficiency losses associated 
with “snowplow” loading of distributed neutral pro- 
pellant. The two-stage scheme used in PPT-1 allows 
the propellant to be loaded in a manner which more 
closely approximates the optimal “slug” loading. 

2. Electrical efficiency: PRC-PPT1 will use a toroidal 
"transmission line” capacitor. The custom-designed 
capacitor, similar to those developed by Hayworth 
et al.[ 14], will be impedance matched to the load 
(the thruster arc discharge) to yield a non-ringing 
current waveform, leading to more efficient trans- 
fer of the capacitively stored energy. The capaci- 
tor design will also lead to a smaller, lighter over- 
all thruster design. The thruster discharge chamber 
will occupy the hollow central region of the capaci- 
tor; the use of a single, large capacitor of this design 
has been shown to be considerably lighter (one-sixth) 
than a conventional parallel combination of equal to- 
tal capacitance! 15]. Also, the thruster geometry has 
been chosen to maximize electrical efficiency. The 
electrode configuration is designed to form a “vir- 
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tual cathode’' (see section 2.1), which will maximize 
the discharge chamber axial inductance gradient and, 
hence, maximize the electrical efficiency [2]. 

3. Electrode erosion: PRC-PPT1 will use a hollow, 
truncated cathode (inner electrode). The aim is to 
distribute the cathode current over a larger area, pre- 
vent current constriction inside the electrode material 
(which occurs at the leading edge of a propagating 
current sheet), and to provide a plasma-electrode in- 
terface that is not starved of propellant. Lithium pro- 
pellant was chosen, for among other reasons, to test 
whether or not the reduced electrode erosion found in 
the Lithium Lorentz Force Accelerator (LiLFA)[16] 
could also be realized in a pulsed plasma thruster. 

4. Reliability: The use of the molten lithium plasma in- 
jector eliminates the need for both a gas valve and a 
high-current electrical switch (both the first and sec- 
ond stage are self-switching); the injector design ful- 
fills both roles using no moving parts to provide, in 
principle, a highly reliable propellant feed and elec- 
trical switching system. The transmission line capac- 
itor eliminates lifetime-lirmting current reversal and 
includes an integrated cooling system to extend ca- 
pacitor life. 

The sections which follow more completely describe 
the geometry and operation of the proposed thruster, and 
present the preliminary design of a proof-of-principle lab- 
oratory thruster. 

2 Thruster Description 

2.1 General description 

The PRC-PPT1 is a two-stage pulsed plasma thruster (the 
use of multi-stage electromagnetic acceleration was first 
advocated by Turchi et ai [17], who used a Teflon-ablative 
PPT to inject plasma into a second, electromagnetic accel- 
erator stage). The components of PRC-PPT1 are schemat- 
ically illustrated in Fig. 1. 

The sequence of events in the operation of the thruster 
is as follows. Liquid lithium is delivered from a reser- 
voir to the first stage using an MHD flow coupler. A 
lithium droplet forms at the outlet of the propellant feed 
system, inside the first stage plasma injector. A high cur- 
rent discharge initiates inside the injector, which vapor- 
izes and ionizes the lithium droplet. The injector ejects 
a dense, thermal plasma into the main discharge cavity 
(second stage). A Townsend avalanche ensues, forming 
a current sheet with mega-amp level current. The current 


sheet accelerates the length of the second stage, entraining 
the propellant and ejecting it from the thruster at a speed 
on the order of 100 km/s. 

Let us consider each of these phases of operation in 
more detail. The propellant feed system is diagrammat- 
ically illustrated in Fig. 2a. Liquid lithium propellant 
is stored in an un- pressurized, heated (~ 1 80 C) reser- 
voir. The propellant is delivered from the reservoir to the 
thruster using an MHD flow coupler[18]. The flow cou- 
pler applies a JxB force to the molten lithium by pass- 
ing current through it, transverse to an applied magnetic 
field. Precise control of the lithium mass flow rate is ac- 
complished by simply adjusting the current that passes 
through the flow coupler. The flow coupler has no mov- 
ing parts - emphasizing an attractive feature of using a 
conductive propellant. 

The lithium is fed into the first stage injector through 
a (electrically insulating) ceramic tube. The first stage 
capacitor is connected between the metal (stainless steel) 
lithium feed line and the outer electrode (anode) of the 
injector. As propellant is fed into the injector, a lithium 
droplet forms at the end of the ceramic feed line, as illus- 
trated in Fig. 2b. As the droplet grows in size, it eventually 
contacts the outer electrode and ”arcs-over” - closing the 
circuit connected to the first stage capacitor. The capaci- 
tor has sufficient energy to vaporize and (singly) ionize the 
entire lithium droplet. The newly formed lithium plasma 
thermally expands out of the first stage injector and into 
the second stage discharge chamber. The vaporization of 
the droplet re-opens the first stage capacitor circuit, al- 
lowing it to be recharged. As lithium is continuously be- 
ing pumped by the flow coupler, a new droplet is formed 
and the process repeats. Again, two more advantages are 
gained through the use of a conductive propellant: no high 
current switch is needed in the first stage, as the circuit 
is self-switching, and the propellant in the feed line acts 
as the cathode, eliminating the first stage cathode erosion 
issue. We are, in essence, feeding a liquid cathode into 
the injector, and using the cathode erosion products as the 
propellant for the second stage. 

The final phase of operation in a single pulse of the 
PRC-PPT1 is the second stage electromagnetic acceler- 
ation. The conceptual (and as yet theoretically and ex- 
perimentally unfounded) evolution of the second stage 
current sheet is envisaged in Fig. 3. The second stage 
electrode configuration, which is essentially a hybrid z- 
pinch/coaxial geometry, uses a truncated inner electrode 
(cathode) to induce a ’’virtual cathode” along the center- 
line of the accelerator, in an attempt to maximize the ac- 
celeration channel inductance gradient and, hence, max- 
imize the force on the current sheet. The current sheet 
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Figure 1 : Schematic ol the PRC-PPT1 geometry and components (sectioned to show inner detail.) 


is expected to initiate near the back of the accelerator - 
where the Poynting vector flux is first able to deposit its 
energy. The J x B force density w ill drive the current sheet 
around the end of the cathode, pinching the plasma at 
the center, in a manner similar to a dense-plasma-focus 
device[?]. Provided that current continues to be driven 
by the capacitor, the anode current sheet attachment will 
continue to propagate axially, and the cathode column will 
be extruded, forming a ’’virtual” center electrode (see Fig. 
3). The cathode current sheet attachment will distribute 
itself over the inside of the hollow cathode, conducting 
through residual lithium plasma from the first stage injec- 
tor. The current sheet will continue to accelerate axially 
until the current subsides, either through exhaustion of the 
capacitively stored energy, or thrc >ugh the disruption of the 
virtual cathode column, from instabilities that are known 
to evolve in z-pinch current configurations. The entrained 
propellant will ultimately be ejected from the thruster, and 
thrust will be derived. 

Other features of the PRC-PPT1 design are illustrated 
in Fig. 1. A metallic UV baffle is located in front of the 
breech insulator. The purpose of the baffle is to shield 
the insulator material from damaging ultraviolet radiation 
which is emitted from the current sheet plasma. Also 


shown are anode and cathode cooling loops. Both elec- 
trodes will be actively cooled using flowing liquid coolant. 

The discussion above describes the general geometry 
and operation of the proposed thruster. The next sec- 
tion gives more specific design details for a first, proof- 
of-principle, thruster. 

2.2 PRC-PPTIa design 

The aim of our first experimental prototype thruster (PRC- 
PPTla) is to validate the operation of the liquid lithium 
plasma feed system, to characterize the ejected plasma 
and test the repetitive operation of the first stage plasma 
injector, and to map out the current sheet trajectory inside 
the second stage of the thruster, to see if the current sheet 
evolution envisaged in Fig. 3 can be realized in an actual 
thruster. 

The design of the prototype thruster is constrained by 
several factors: first, of course, the thruster is expected to 
meet the performance goals stated in section 1.2, second, 
the components must be constructed from materials on- 
hand at the PRC and, lastly, the operation of the thruster 
must be compatible with vacuum facilities presently avail- 
able at the PRC. The last two constraints limit us to, at the 
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a) first stage 



MHD flow coupler 


b) 


© droplet formation 



© droplet- anode contact 



© arc discharge and plasma ejection 


1. initiation 



2. pinching 



3. axial acceleration 




Figure 2: PRC-PPT1 propellant teed system: a) diagram 
of liquid lithium feed system, b) first stage droplet forma- 
tion and plasma ejection sequence. 


Figure 3: Conceptual evolution of current sheet in second 
stage of PRC-PPT1: initiation, pinching, and acceleration 
with virtual cathode. 


present time, the design of a thruster that operates in a 
single shot mode, as we presently do not posses: capac- 
itors suitable for repetitive operation, a high-power high- 
voltage power supply, or a vacuum facility capable of han- 
dling large mass flow rates. Components that are readily 
available in our laboratory inventory include 100 17.5 /xF, 
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40 kV, low inductance (30 nH) capacitors, which were re- 
cently acquired from General Atomics Inc.. The capaci- 
tors have single ended coaxial connectors to facilitate low 
inductance connection to the thruster. A large selection of 
raw copper (including 20 cm diameter thick- walled pipe) 
material is available for use as electrode material. 

A schematic illustration of the PRC-PPTla complete 
assembly and a cross-section of the discharge chamber is 
shown in Fig. 4. The thruster will use four of the afore- 
mentioned capacitors, connected in parallel, to allow up 
to 50 kJ discharge energy. A cross-sectional schematic of 
the discharge chamber is shown in Fig. 4b, where the rel- 
evant dimensions are shown. In the illustration, only the 
anode inner radius is explicitly given; this is the dimen- 
sion of the copper pipe that is available to us. The other 
dimensions, L a , L c , and R c (the anode length, the cath- 
ode length, and the cathode radius, respectively) are to be 
established through modelling. 

Since the dimensions of the thruster will be established 
through preliminary modelling efforts, let us list, in one 
place, the type of modelling efforts that will be needed to 
guide the physical dimensioning of the prototype thruster: 

1. Lithium droplet decomposition. Modelling of the 
lithium droplet vaporization in the first stage plasma 
injector is necessary to predict the temporal evolu- 
tion and thermodynamic slate of the decomposition 
products, as well as to determine the first stage ca- 
pacitor energy requirement. A model of the droplet 
formation itself will be necessary to size the first 
stage liquid injection orifice and electrode spacing. 

2. First stage injection. The expansion of the first 
stage injector plasma should be modelled to guide 
the design of the injector and the second stage cath- 
ode. It is imperative that the first stage injector 
and second stage cathode deliver the majority of 
the lithium plasma propellant into the second stage 
discharge chamber before the second stage arc dis- 
charge commences, otherwise propellant will get left 
behind the accelerating current sheet and, hence, 
constitute a propellant utilisation inefficiency. A 2-d 
fluid model would be helpful in estimating the tem- 
poral evolution of the lithium plasma density profile 
as it thermally expands out of the hollow cathode. A 
kinetic model will be required to accurately estimate 
the the time of breakdown. 

3. Current sheet propagation. Modelling of the cur- 
rent sheet propagation is necessary to determine the 
proper length of the second stage anode. A 0-d cur- 
rent sheet/lumped circuit element model would pro- 
vide a first estimate. A full 2-d MHD simulation 
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would provide more accurate results, as well as elu- 
cidate the expected current sheet configuration, that 
is, show how the current sheet tilts as it propagates. 
This information could be used to tailor the first stage 
plasma injector to provide an initial mass distribution 
which inhibits adverse tilting of the current sheet. 

4. Stability of the virtual cathode. The PRC-PPT1 

design uses a hybrid z-pinch/coaxial electrode con- 
figuration, which is intended to produce a propa- 
gating planar current sheet electrically connected to 
the (physical) cathode at the back of the accelera- 
tor through a long plasma column (the virtual cath- 
ode). The virtual cathode is essentially a z-pinch cur- 
rent configuration, which is known to be unstable. 
For the purposes of the PRC-PPT1 design, a mod- 
elling effort, which draws upon existing knowledge, 
is needed to predict how long the virtual cathode can 
be maintained inside the accelerator, before current 
flow is disrupted by current driven instabilities. The 
onset of such instabilities may limit the length of the 
acceleration channel that can be practically be imple- 
mented. 

As stated in the introduction, the present work is pri- 
marily conceptual in nature - the necessary modelling 
effort enumerated above has not been completed. To 
date, the first three items enumerated above have been ad- 
dressed to a limited extent; the results of those efforts are 
presented in the next section. 

3 Modelling 

3.1 Current sheet propagation 

A 0-d plasma thruster code (PTC1) was developed to ob- 
tain a rough estimate of the expected current sheet tra- 
jectory. The code is used to solve the unsteady thruster 
electrical circuit equations, coupled with a 1-d model of 
the current sheet dynamics. A schematic of the compu- 
tational elements is shown in Fig. 5. The lumped circuit 
elements represent the following plasma thruster compo- 
nents: 

C = capacitor capacitance 

R1 = transmission line resistance 

R2 = “crowbar” resistance 

R3 = current sheet resistance 

LI = capacitor internal inductance 

L2 = transmission line inductance 

L3 = instantaneous thruster inductance 
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a) 


b) 


capacitors (4) 


+ 





Figure 4: a) Isometric view of complete PRC-PPTla assembly, b) Cross-section of PRC-PPT la discharge chamber. 


Applying Kirchoff’s circuit laws and Newton’s second 
law to the circuit illustrated in Fig. 5 results in five first- 
order simultaneous ODEs: 



Q 

= -/l . 


(1) 

1 

L\ T L/2 

Q 

<7 

(R\ + #2 )h 

+ 

, (2) 

h = ~ [R 2 I 1 - 

^3 

- (R 2 4 R3 T 

L'v)h ] 

, O) 
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(4) 


V — 

1 r'r 2 
2m 


(5) 


where Q is the instantaneous charge on the capacitor, 
V is the inductance gradient in the gun, v is the cur- 
rent sheet speed, x is the instantaneous current sheet ax- 
ial position, and m is the instantaneous mass in the cur- 
rent sheet. Coaxial electrodes are assumed with constant 



R1 L2 




Figure 5: Computational circuit elements of PTC 1. 


cross-section (no taper), in which case the inductance gra- 
dient is constant[19] 

In — , (6) 

2n Ti 

where fi 0 is the permeability of free space, and r Q and r* 
are the outer and inner electrode radii, respectively. The 
instantaneous current sheet mass is calculated using either 
a slug or snowplow model, as specified by the user. Equa- 
tions 1-5 were integrated using Matlab. 

The following circuit parameter values were used to 
simulate the performance of PRC-PPTla: 

C = 70 fiF 

ri = o n 

R2= 100 ft 
R3 = 2.5 x 10” 3 ft 
LI =7.5 x 10" 9 H 
L2 = 20 x 10" 9 H 

An initial voltage of 37.8 kV was prescribed (E = 50 kJ) 
with 5 mg of lithium propellant loaded as a slug. The 
outer electrode radius was set to 10 cm and the inner elec- 
trode radius was set to 1 cm (the virtual cathode is as- 
sumed to form a thin filament of current along the center- 
line; a typical current sheet thickness is assumed[13]). 

Figure 6 shows the results of the simulation. The cur- 
rent waveform is predicted to be over-damped, reaching 
a peak of about 1 MA with a pulse length of about 6 /xs. 
The current sheet is predicted to reach a terminal speed 
of about 1 10 km/s over about 50 cm of propagation. The 
total impulse bit derived during a single pulse is estimated 
to be about 0.6 N-s. 

The results of the 0-d model indicate that PRC-PPTla 
will use on the order of 5 mg of lithium per pulse, and re- 
quire an anode length (L a ) of about 50 cm. The oversim- 
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Figure 6: Results of 0-d analysis: current, sheet speed, 
sheet position, and thrust as a function of time. 

plifications inherent in the 0-d model will almost certainly 
lead to an over-prediction of performance; nevertheless, 
the results do give us a starting point to begin more so- 


phisticated analysis. 

3.2 Lithium droplet decomposition 

Based on the results of the 0-d model we can begin to 
estimate the size of the first stage propellant feed system 
and capacitor energy requirement. 

The density of liquid lithium is 515 kg m -3 [20]. Thus, 
assuming a 5 mg/droplet, the diameter of the lithium 
droplet which must be introduced into the injector is 2.6 
mm. 

The energy required to vaporize and (singly)ionize the 
lithium droplet may be estimated by calculating the num- 
ber of moles of lithium atoms in a 5 mg sample and multi- 
plying by the heat of vaporization (134.7 kJ/mol) and the 
ionization energy (513.3 kJ/mol)[20]. The calculation in- 
dicates that about 460 J of energy is needed, sizing the 
first stage capacitor energy to at least 0.5 kJ. 

3.3 First stage injection 

The injection of the lithium plasma from the first stage 
into the second stage discharge chamber was modelled 
using the 2-D MHD code MACH2[21]. A computational 
mesh which uses the PRC-PPTla electrode geometry was 
implemented. The thermal expansion of 5 mg of lithium 
propellant (1 eV initial temperature) from the inside of the 
hollow cathode was simulated. The computational mesh 
and the spatial evolution of the injected plasma at three 
different times is illustrated in Fig. 7. 

The results of the simulation indicate that about ninety 
percent of the propellant is evacuated from inside the hol- 
low cathode in about 10 ps. NOTE: more material will be 
added here . 

The use of such a powerful and complicated numerical 
tool as MACH2 to simulate the first stage plasma injec- 
tion may seem to be somewhat of an “overkill”. How- 
ever, we are presently in the process of implementing the 
electromagnetic and the circuit simulating capabilities of 
MACH2; the propellant loading illustrated in Fig. 7 will 
serve as an initial condition for our full, unsteady MHD 
simulations. 

4 Conclusion 

Design details of a proposed high-energy (~50 kJ/pulse), 
two-stage pulsed plasma thruster were presented. The 
long-term goal of this project is to develop a high-power 
(~500 kW), high specific impulse (~7500 s), highly effi- 
cient (~50%), and mechanically simple thruster for use 
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Figure 7: Spatial evolution of the injected plasma at 
three different times (figure needs some additional work 
-TEM. 

as primary propulsion in a high -power nuclear electric 
propulsion system. The proposed thruster (PRC-PPT1) 
utilizes a valveless, liquid lithium-fed thermal plasma in- 
jector (first stage) followed by a high-energy pulsed elec- 
tromagnetic accelerator (second stage). A numerical cir- 
cuit model coupled with one-dimensional current sheet 
dynamics, as well as a two-dimensional numerical MHD 
simulation, was used to qualitatively predict the thermal 
plasma injection and current sheet dynamics, as well as to 
estimate the projected performance of the thruster. The re- 
sults of the 0-d circuit model indicate that five milligrams 
of propellant can be accelerated to 1 100s in a half-meter 
long accelerator, using a capacitor bank energy of 50 kJ. 
While the 0-d code is likely to grossly over-predict the 
performance of an actual thruster, it provides a starting 
point to begin more sophisticated analysis. 


The use of a molten metal propellant and a two stage 
acceleration scheme is shown to provide potential benefits 
which may be exploited to overcome many of the negative 
issues associated with PPTs, such as propellant utilization 
inefficiency and mechanical reliability. A set of further 
modelling efforts, and the experimental testing of a proto- 
type thruster, is suggested to determine the feasibility of 
demonstrating a full scale high-power thruster. 
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Why do we need electromagnetic accelerators? 



■ Electrothermal acceleration 

■ Low specific impulse, high thrust efficiency, high thrust density 

■ Electrostatic acceleration 

■ High specific impulse, high thrust efficiency, low thrust density 

■ Electromagnetic acceleration 

■ Steady-flow accelerators (Magnetoplasmadynamic (MPD) thruster) 

■ Intermediate specific impulse, high thrust efficiency, high thrust density 

■ Pulsed accelerators (Pulsed Plasma Thruster (PPT), QS-MPD) 

■ High specific impulse, low thrust efficiency, high thrust density 
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(From R.G. Jahn, Physics of Electric Propulsion) 
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Quasi-steady MPD versus Pulsed Plasma Thruster 




■ PPTs offer several potential advantages over Quasi- 
steady MPDs 

■ Higher thrust efficiency due to higher current density 

■ Higher average thrust density 

■ Lower electrode erosion 

■ Reduced frozen flow losses 

■ More effective coupling between accelerating field and ambient 
(propellant) gas 

■ Lower system mass (single capacitor versus a PFN) 


feJHfr£A-;.lvy 
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PPT Problems and Solutions 




Propulsion Rososrch Center 


Problem 

Solution 

■ Efficiency 

■ Efficiency 

■ Poor propellant utilization 
efficiencv 

■ Liquid propellant injection using DENSE 
PLASMA INJECTORS 

■ Poor dynamic efficiency 

■ Slug-like propellant loading using DENSE 
PLASMA INJECTORS 

■ Poor electrical efficiency 

■ Impedance matching using 

TRANSMISSION LINE CAPACITOR 

■ Electrode Erosion 

■ Electrode Erosion 

■ Limited data 

■ COMPREHENSIVE EROSION STUDY 

■ Cathode erosion 

■ High density plasma using DENSE 
PLASMA INJECTORS 


■ NOVEL THRUSTER GEOMETRY 

■ Reliability 

■ Reliability 

■ High current switch 

■ Eliminate switch using DENSE PLASMA 
INJECTORS 

■ High speed gas valve 

■ Liquid propellant injection using DENSE 
PLASMA INJECTORS 
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HPPPTX PPT-1 (schematic) 



Propulsion ftosovcti Contv 


Target Performance 
Isp ~ 9000 s 
r) t > 50 % 

T/A ~ 10 5 -10 6 N/m 2 


"transmission line" capacitor 



US 
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> 1 st Stage: Theta pinch 

> E vap = 0.02 J 

> E ion = 74.7 J 

> E cap ~ 750 J 
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DESIGN CALCULATIONS 
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HPPPTX Program Goals 




■ Phase I 

■ Develop DENSE PLASMA INJECTOR 

■ Develop TRANSMISSION LINE CAPACITOR 

■ Develop HPPPTX PPT-1 thruster 

■ Phase II 

■ Operate and measure performance of HPPPTX PPT-1 at 100 kW 
CW power level 

■ Conduct comprehensive electrode erosion study 

■ Develop HPPPTX PPT-2 thruster 

■ Phase III 

■ Operate and measure performance of HPPPTX PPT-2 at 1 MW 
CW power level (world record) 

■ Advance TRL level for High Power Plasma 
Propulsion 
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